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1. 


Statement  of  Work 


Engineering  analysis  must  be  performed  to  investigate  the  implementation  of  a  sliding 
mode  controller  algorithm  that  will  reduce  sensitivity  to  feedback  noise  in  Control  Actuation 
Systems  (CAS)  such  as  those  used  on  Guided  MLRS,  Low  Cost  Precision  Kill  (LCPK),  and 
Compact  Kinetic  Energy  Missile  (CKEM).  The  GEE  model  of  the  CAS  shall  be  used  to 
demonstrate  performance  enhancement  achieved  with  the  sliding  mode  controller  for  a  variety  of 
potential  noise  sources. 
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2.  Plant  Description  and  Problem  Formulation 


The  control  problem,  which  is  addressed  in  this  technical  report,  consists  of  controlling  a 
nonlinear  servomechanism  affected  by  noise.  A  block-diagram  of  the  servomechanism  with 
disturbance  and  noise  are  shown  in  figures  1  and  2 


Fig.  1  A  block-diagram  of  the  servomechanism 
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j  Annatme  vxjitage  (Volt) 


dt 

Denvaiivr 


jKQise(Volt)| 


Fig.  2  A  block-diagram  of  disturbance  and  noise 


Neglecting  backlash  and  hysteresis  on  the  step  of  a  sliding  mode  controller  design  the 
state  variable  model  of  the  servomechanism  (figures  1  and  2)  is  derived  as  follows: 

1 

— -  =  — 0}„ 
dt  n, 

,  (1) 

+  "v  -  RJa  -  ] 


where 
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0^  is  position  of  the  load  {rad),  is  an  armature  current  (A),  ci),„  is  motor  angular  velocity 
(rad/s) ,  Hg  =  125.9  an  actuator  gear  ratio.  6^  =  is  position  of  the  motor  (rad),  |v^|  <  31  is 
armature  voltage  (V),  is  voltage  noise  (V),  k,  =  0.445  is  a  torque  constant  (in-lbflA) ,  is 
back-EMF  constant  (V/rad/s),  R^=\0Q.  is  motor  winding  resistance.  =0.0036^  is  motor 
winding  inductance,  7,  =0.0000125  is  total  inertia  reflected  to  the  motor  {in-lbf  -s^). 


7^=0.003  is  load  inertia  {in-lbf  -s'),  7  =  7, =  0.0000127  is  total  inertia 

n: 


(in  -Ibf-s^),  F{.)  =  32ie^  +  0. 1 8  signd„^  [0.8 1 8|0^  |  +  (o.8 1 8|0^  |  + 1)'  f  is  disturbance  torque 


{in-lbf). 

The  servomechanism  control  problem  is  to  design  a  control  law  in  terms  of  the  armature 
voltage  for  the  servomechanism  to  achieve  accurate  tracking  of  the  input  command  profile  diJt) 

in  terms  of  load  position,  which  is  generated  by  the  guidance  system  in  a  real  time,  regardless  the 
voltage  noise  and  the  torque  disturbance  F{.) . 


3.  Basics  of  SISO  Nonlinear  Output  Tracking  on  Sliding  Modes  in  Canonical  Space 

A  nonlinear  MIMO  system  s  considered 

\x  =  f(x,t)  +  A/fJC,  t)  +  [B(x,  t)  +  AB(x,  t)]u 
[y  =  h(x) 

where  xe  R".  ye  R\  ueR\  f(x,t),  Af{x,t),  B(x,t)  and  AB(x,t)  are  smooth  functions.  Here 
Af(x,t)  and  AB(x,t)  represent  system  uncertainties  and  disturbances.  A  nominal  system  is 
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obtained  when  A/T-i;tj=0  and  Afi('x,r)  =  0.  Tlie  problem  is  to  design  a  control  function  u  in 
order  to  achieve  a  robust  asymptotic  tracking  the  command  output  profile  yjt)  V/  =  1^: 

\\me  =  Vm[y^(t)-y(t)]=0.  (3) 

A  sliding  mode  controller  (SMC)  design  consists  of  two  steps  [1-7]. 

In  the  first  step,  a  sliding  surface  is  chosen  such  that  the  closed  loop  system  motion  on  this 
surface  exhibits  desired  behavior  regardless  of  plant  uncertainties,  and  disturbances. 

In  the  second  step,  a  control  function  is  chosen  to  provide  reaching  of  the  sliding  surface  by  the 
system  state  in  finite  time  and  to  guarantee  system  motion  in  this  surface  thereafter.  System 
motion  on  the  sliding  surface  is  called  a  sliding  mode,  and  strict  enforcement  of  the  sliding  mode 
typically  leads  to  discontinuous  high  frequency  switching  control  functions. 

3.1  System  transformation 

Assuming  the  system's  (2)  relative  degree  equal  to  the  system's  order:  r  =  n,  the  system  (2)  can 
be  easily  transformed  to  a  companion  form  (or  canonical  space)  [7]  as  follows: 

^2’  ^2  f)  +  Aipf.Y. r) -(■  r)  +  Afefx, r)]M,  y  =  ^j,  (4) 

where  (p(x,t)+A^x,t)  =  L"f(h),  b(x.t)  +  Ab(x,t)=LJL”f^(h))  and  L"f(.),LJ.)  are  corresponding 
Lie  derivatives.  It  is  assumed  that  bix,t)  +  Abix,t)^0  in  a  reasonable  domain  xsF.  The 
functions  A(p(x,t)  and  Ab(x,t)  reflect  uncertainties  and  disturbances  in  the  original  system  (2). 

3.2  Sliding  surface  design 

The  switching  surface  can  be  designed  [1-7]  as  follows: 
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(5) 


(J  =  ® ' 

(=0 


where  r=n  is  a  relative  degree  of  the  system  (2).  The  coefficients  c,.V/  =  l,n-l  must  be 

identified  to  provide  the  desired  eigenvalues  placement  to  the  differential  equation  (5).  It  is 
obvious  that  the  systems  motion  on  the  switching  surface  (5)  implies  the  robust  output  tracking. 
3.3  Control  function  design 

A  Lyapunov-based  design  approach  is  used  to  derive  a  control  law.  The  goal  is  to  design 
such  control  function  u  that  the  sliding  surface  (5)  will  be  reached  in  a  finite  time.  Assuming 
c„_i  =1  the  CT  -  dynamics  is  derived  as  follows: 

o  =  cr(0)  =  cro,  (6) 

where  <D(.)  =  +...  +  Coe-(p(x,t)-A(p(x,t). 

Global  asymptotic  stability  of  the  equilibrium  point  cr  =  0  of  the  equation  (6)  will  be  achieved 
by  a  control  law,  which  is  designed  using  Lyapunov  function  technique.  Assuming  b(x,t)  >  0  in  a 


reasonable  domain,  the  following  candidate  to  a  Lyapunov  function  is  used 


(7) 


and  its  derivative  is  calculated  on  the  trajectory  of  the  system  (6)  as  follows: 

Q  =  QbQcr  +  b-^  (.)d>(.)  -  (1  +  5  )u  j ,  (8) 

.  -  A*(.) 

where  o  —  .  In  order  to  provide  asymptotic  stability  to  the  origin  of  the  system  (6)  the 

U\*) 

following  format  for  the  derivative  of  the  candidate  to  the  Lyapunov  function  is  required 
2<-PiK  Pi>0-  (9) 
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It  is  well  known  that  if  inequality  (9)  holds  then  the  sliding  surface  (5)  will  be  reached  in  a  finite 
time,  which  can  be  calculated  as  follows: 


|g(0)| 

Pi 


(10) 


A  following  control  law  is  designed  to  asymptotically  stabilize  cr  =  0  in  equation  (6) 

u=u^^  +  p signer,  p>0,  (11) 

where  is  an  estimate  the  equivalent  control  function  that  makes  derivative  of  the 
candidate  to  the  Lyapunov  function  equal  to  zero.  The  equivalent  control  function  is 
identified  from  equation  (8).  This  is 


(12) 


Remark  1.  We  have  to  use  instead  of  in  the  control  law  (11).  The  reason  is  that  the 

equivalent  control  (12)  cannot  be  implemented  because  it  contains  uncertain  terms  Ab{.) ,  A<p{.) . 

Assume  that  the  estimattion  eiTor  of  the  equivalent  control  is  bounded  in  a  reasonable 
domain.  This  is 

(13) 


1,- 


where  Am,,  =  i.)bi-)<r  +  i’"*  (.)<!>(•)  -  (l  +  5 )m„  . 


Substituting  equation  (11)  into  equation  (8)  we  can  obtain  the  following  expression: 

0  =  cr[AM„  -  (l  +  5 )psign(7]  <  |tT|[M  -  (1  +  5)p] .  (14) 

Obviously  inequality  (14)  will  fit  inequality  (9)  if 
M  +  p, 


P> 


1  +  5 


(15) 
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Remark  2.  Inequality  (15)  will  lead  to  a  positive  p  if  ISl  =  <  1 .  Tliis  condition  limits  the 

'  '  hi) 

value  of  uncertainty  which  the  sliding  mode  controller  (5),  (11)  and  (15)  can  tolerate. 

Remark  3.  The  sliding  mode  controller  (5),  (11)  and  (15)  is  robust,  because  the  control 


amplitudes  are  selected  on  the  basis  of  inequalities.  In  this  case  it  is  sufficient  to  know  only 


ranges  of  elements  of  uncertain  functions  <  M  and 


Remark  4.  The  downside  of  the  proposed  approach  is  in  necessity  of  calculation  (estimation)  of 
canonical  variables  of  the  system  (4),  It  practically  means  that  we  need  to  measure  or  estimate 


controlled  output  and  its  derivatives.  It  can  be  done  using  asymptotic  observers  or  sliding  mode 
observers. 


3.4  Chattering 

In  order  to  avoid  a  control  chattering  a  continuous  implementation  of  a  discontinuous 
control  law  is  often  considered  [7-9].  Usually  this  is  an  equivalent  control  plus  a  linear  control 
law  with  respect  to  sliding  surface  or  state  variables.  Very  often  a  linear  control  law  is  taken  with 

saturation  sat  — ,  where  £,.  is  a  thickness  of  a  boundaiy  layer  in  i  th  switching  element  [7-9]. 

k  ) 

Such  sliding  mode  controllers  reduce  the  control  activity,  provide  a  smoother  operation  and  are 
free  from  chattering,  but  have  some  disadvantages  comparing  to  infinity  frequency  switching 
sliding  mode  controllers.  They  are  [7-9]: 

1.  Sliding  surface  can  be  reached  only  asymptotically,  i.e.  reaching  time  «= . 
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2.  In  case  of  nonvanishing  disturbances  and  uncertainties  convergence  is  provided  only  to  £,  - 

vicinity  of  the  sliding  surface.  A  finite  reaching  time  to  this  vicinity  is  guaranteed.. 

Commenting  the  first  disadvantage  we  can  say  that  instead  of  reaching  time  we  can 
operate  with  settling  time  [8,9],  which  is  finite  and  can  be  controlled  by  a  con'esponding  choice 
of  parameters  of  the  sliding  surfaces  (11).  Commenting  the  second  disadvantage  we  can  say  that 
the  distance  between  the  system's  trajectory  belonging  to  the  e,  -  vicinity  of  the  sliding  surface 

and  the  system's  trajectory  belonging  to  the  sliding  surfaces  proportional  to  the  size  of  this 
vicinity  [1,2].  In  this  case  we  can  talk  about  e,  -  robustness  of  the  sliding  mode  with  a 

continuous  control  signal.  The  price  to  pay  is  an  eventual  steady  state  eiTor  while  the  controller 
operates  lineaily  in  the  e,  -  vicinity  of  the  sliding  surface.  This  problem  can  be  addressed  by  the 

introduction  of  an  integral  term  into  the  equation  of  the  sliding  surface  (11).  The  continuous 
implementations  of  a  sliding  mode  controller  are  widely  known,  for  instance  [10-15]. 

4.  Basics  of  nonlinear  output  tracking  on  sliding  modes  in  physical  state  space  [4] 

It  is  well  known  that  control  systems  with  sliding  modes  (SM)  are  robust  to  parameter 
uncertainties  and  disturbances  which  makes  them  attractive  for  output  tracking  in  MIMO  systems 
[1-3].  In  order  to  achieve  the  desired  output  tracking,  high  order  derivatives  of  the  output 
tracking  error  are  generally  required  [7].  In  order  to  eliminate  output  error,  multiple 
differentiation  coordinate  transformations  were  developed  [7].  The  basis  transformation  proposed 
in  this  work  identifies  such  “prescribed”  profiles  for  new  state  (physical)  variables  so  that  the 
state  tracking  yields  the  output  tracking,  while  unmatched  disturbances  and  nonlinearities  of  the 
original  system  become  matched  in  the  new  basis.  Conventional  and  dynamic  sliding  mode 
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controllers  (SMC)  are  designed  to  provide  the  desired  output  tracking  in  a  SM.  The  contribution 
of  this  work  consists  in 

•  the  development  of  a  nonlinear  transfoitnation  of  the  original  system  with  unmatched 
disturbances  and  nonlinearities  into  a  form  convenient  for  a  SMC  design, 

•  conventional  and  dynamic  SMCs  design  for  the  desired  output  tracking.  Both  conventional 
and  dynamic  SMCs  are  insensitive  to  matched  disturbances,  but  the  dynamic  SMC 
accommodates  to  unmatched  ones. 

4.1  Problem  formulation 

The  following  plant  is  considered: 

x=  Ax  +  F(x,t}-^Bu,  y  =  Gx,  (16) 

where  xeR"  -  a.  state  vector;  ueR'"  -  a  control  function;  ysR”'  -  an  output;  A,B,G  -  constant 
matrices;  F(x,t)<aR''  -  an  unmatched  disturbance.  {A,B}-  a  controllable  pair.  rankB  =  m, 
y’  it)  e  /?”-  reference  output  profile.  We  wish  to  specify  a  SMC 
\ut,<5>.(x,y\t)>0, 

U;  =  <  (17) 

[u:.  O;  (x,  y\ t)  <  0,  Vi  =  I,  m, 

where  (x,y  ,t)  are  conventional  sliding  manifolds  or  sliding  manifolds  as  dynamic 

operators,  representable  as  transfer  functions  acting  on  some  states  and  on  output  tracking  errors 
of  the  system  (1);  i/*, -  continuous  functions  of  x  and  t ,  to  accomplish  the  following  goals: 

•  The  motion  of  the  output  tracking  eiror  e(t)  =  y’(t)-y(t)  in  the  SM  is  described  by 
homogeneous  time-invariant  linear  differential  equations  with  given  eigenvalues  placement. 

•  The  existence  of  the  SM  in  the  designed  sliding  manifolds  must  be  provided. 
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4.2.  Nonlinear  transformation  of  the  system  basis 


In  order  to  eliminate  the  output  enor  multiple  differentiation  [7],  the  output  tracking  is 
transformed  to  an  equivalent  state  control  problem.  Tlie  following  lemma  is  proved. 

Lemma.  The  system  (1)  in  the  new  basis 

(18) 

has  the  following  convenient  form: 

=A^‘z‘  +A“z^  -B^u  +[x-'  -A^‘x‘'  -A^\x-'  -FUx,t)] 
e  =  G‘z‘+G^z-:  : 


b'l 

1 

# 

_ 1 

1 - 

1 _ 

1 

1 

« 

- 1 

(19) 


R"-'":  z^eeR"-, 


where 


nonsingular  linear  transformation  M  = 


'm^ 

M' 


has  been  specified  as  follows  [1,2]: 


MB  = 


MAM-‘  = 


'm^b' 

'0  ' 

M^B 

L  det  B  ^  0 . 


'a"  A^' 
A^^A^^ 


'M‘ 


F  = 


F' 

F^ 


'm‘ 


-1-/ 


=  [g‘\g']. 


(20) 


(21) 


x‘*  ,x^'  satisfy  the  nonlinear  system  of  differential-algebraic  equations  (DAE)  [4]: 


x^* 

"a'' 

'x‘*' 

’fUxJ)' 

1 - 

o 

1 _ 

G‘ 

G\ 

y\ 

_-y’(0 . 

(22) 
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Remark  1. 


•  The  reference  state  variable  profiles  x‘*  ,x'*  are  similar  to  “vinual  control  laws”  in  the 
“backstepping”  method  [16],  However,  x‘*  and  .r'*  are  identified  in  one  step  through  a 
solution  of  the  system  of  DAE  (22),  which  is  named  the  system  center  [4]. 

•  Assuming  the  unmatched  disturbance  F'(x,  t)  is  estimated,  the  system  (22)  is  solvable,  if  and 


only  if,  the  matrix  pencil  A 


o' 

-h 

■a“ 

A'"’ 

_0 

0_ 

G^ 

G\ 

.  A  e  is  a  regular  one  [4]. 


•  Implementing  the  concept  of  the  "computed  torque  control,"  the  system  of  DAE  (22)  is 
rewritten  as  follows: 


x‘^-A‘^x‘*-A‘fx^'-F‘(x\t)  =  0,  y'-G^x^*-G^x^'  =0,  x  = 


(23) 


and  only  the  time-dependent  part  of  the  unmatched  disturbance  F'(x*,  t)  must  be  estimated. 


The  transformation  (20)  exists  but  is  not  unique  if  the  pair  {A,  B]  is  controllable  [1,2], 


4.3  Design  of  the  conventional  SMC 

The  following  theorem  is  formulated  and  proved. 

Theorem  1.  Suppose  that  the  SM  exists  in  a  conventional  sliding  manifold 
G  =  z^+cz‘=0,  cr€/?'",  C  e  ,  (24) 

then  the  output  tracking  error  e( t  )  =  y* ( t )- y( t )  of  tht  system  (19)  in  the  sliding  manifold  (24)  is 
described  by  the  following  linear  time-invariant  homogeneous  system  of  differential  equations: 

=  (A“  -  A‘^c)z‘ .  e=(G‘-G\)z‘.  (25) 

Proof.  Substituting  the  equation  (24)  into  the  equations  (19)  yields  the  equations  (25).  ■ 
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Remark  2.  It  is  well  known  [1.2]  that  if  the  pair 


in  the  system  (19)  is 


K.. 

-  i 

■  0  ' 

A^‘A-  _ 

} 

-B\ 

controllable,  then  the  pair  A“,A‘'  in  the  system  (25)  is  also  controllable.  Consequently,  the 
matrix  c  can  be  designed  to  provide  the  given  eigenvalues  placement  to  the  matrix  A“  -  A‘^c. 

The  control  tunction  (17)  must  be  designed  to  provide  the  existence  of  the  SM  in  the 
system  (4)  in  the  conventional  sliding  manifold  (24).  The  following  theorem  is  proved. 

Theorem  2.  The  SM  exists  in  the  system  (19)  in  the  conventional  sliding  manifold  (24)  if 
u  =  u^^+(B^YsIGN{g)-R,  (26) 

where  SIGN(a)  =  diag[sign  =  p,  >  0  i  =  l,m,  and  is  identified  as 

+  cA^')z‘  +(A^^  +  cA^^)z^  +0],  O  =  -  F'(x,t).  (27) 

Proof.  A  candidate  to  the  Lyapunov's  function 

V  =  0.5  g^  a  >  0,  (28) 

was  built,  and  its  derivative  was  identified  as  follows: 

V  =  (z^  +cz')^[(A'^^  +cA^‘)z'  +(A“  +cA^^)z^  (29) 

Substituting  the  formulas  (26),  (27)  into  the  expression  (29)  3aelds  the  following  inequality: 


m 

V  =  -G^  ■SIGN{g)R=  -Y,Pk-Yk\<0  Vp,  >0  \fk  =  T^. 


(30) 


k=l 


The  equilibrium  point  c  =  0  is  asymptotically  stable.  The  SM  exists  in  the  surface  (24).  ■ 

Remark  3.  Analyzing  robustness  of  the  SMC  (24),  (26),  assume  the  equivalent  control  is 


estimated  as  and  ,  then  = 


(=1 


-P,+ 


\  j=i 


signG; 


,  where 
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bij  V/  =  l.n  y/  =  Lin  are  the  components  of  the  matrix  .  and  are  the  components  of  the 


vector  Au^^.  Apparently  the  SM  still  exists  if  p,  >  max 


V/  =  1,/z . 


4.4  Design  of  the  dynamic  SMC 

Assume  the  unmatched  disturbance  F'(x,t)  is  estimated  as  F‘(x,t),  and 
AF‘(x,t)  =  F‘(x,t)-  F‘(x,t).  Then  the  system  (16)  is  described  in  the  new  basis 
'z‘  =A“z'  -AF‘(x,t), 

<e  =A^‘z'  +A^h"  -B^u  +[.e^ -A^‘x'' -A‘\x-' -FUxJ)],  (31) 

e  =  G‘z‘+G^z^\  z^sR"-"':  z^eeR"”. 


In  order  to  reduce  the  effect  of  a  remaining  part  of  the  unmatched  disturbance  AF‘(x,t)  to 


the  output  steady  state  tracking  error,  the  dynamic  sliding  manifold  [4-6]  was  designed  as  a 
linear  dynamic  operator  acting  on  some  states  and  on  output  tracking  errors  [4-6]  as  follows: 


Z(z\e)  =  z^  +  W(s)e  =  [/;  • 


(32) 


where 

s  =  d/dt,  TO  =  [iy,/^)],  W,^(s)  =  P„,Js)/Q„Js)  '^i  =  I,m^j=l.m; 


P„,M  Qn:M) 


are  polynomials  of  s  and  m,  ^ .  The  following  theorem  is  formulated  and  proved. 

Theorem  3.  Suppose  that  the  SM  exists  in  the  system  (31)  in  the  sliding  manifold 
3(z\e)  =  0.  (33) 

then  the  output  tracking  error  e  in  the  system  (31)  in  the  dynamic  sliding  manifold  (33)  is 
described  by  the  operator  equation 
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[/  +  G ' W(s)  +  G ‘ (si  -A“ )-'  A e  =  -G‘ (si  -  A“  r  AF' (x.  t) .  (34) 

Proof.  Assume  the  SM  exists  in  the  system  (31)  in  the  dynamic  sliding  manifold  (33),  and 
substitute  (32),  (33)  into  (31)  to  obtain  (34). 

Remark  4.  In  order  to  provide  asymptotic  output  tracking  in  the  SM.  the  operator  polynomials 
P,n^.U)’Q„..(>^)  =  l,myj  =  I,in  must  be  identified  to  guarantee  the  given  eigenvalues 

placement  and  to  reject  the  effect  of  the  disturbance  AF‘(x,t)  on  the  output  steady  state  tracking 
error.  The  nonlinear  function  AF'(x,t)  can  be  approximated  by  the  time-function  AFfx*(t),t), 
where  .v*(r)  is  the  reference  state  of  (23).  Tlie  disturbance  rejection  condition  lim  sE(s)  =  0  can 

f— >«» 

be  easily  met  by  a  coiresponding  choice  of  the  polynomials  (s),  Q„,  (s) .  The  control  function 

(17)  must  be  designed  to  ensure  the  existence  of  the  SM  in  the  system  (3 1)  in  the  dynamic  sliding 
manifold  (33).  The  following  theorem  is  formulated  and  proved. 

Theorem  4.  The  SM  exists  in  the  system  (31)  in  the  dynamic  sliding  manifold  (33),  if 

u  =  u,^+{b^Y  ■(!+  W(.s)G^)''  •  SIGN{Z) ■  T ,  (35) 

where  SIGN(S)  =  diag{  sign  3;}  and  T  =  [p^,p2,....p„]^.  p,.  >  0  i=  l,m  are  arbitrary  chosen 

constants,  and  is  the  following  equivalent  control  function  [1.2]; 

Uei  =  ■  (7+ WT^)Gj)“'  •  [(A^‘  +  W(s)G‘A“  +  W(s)G^A^^)z‘  +(A‘"  +  W(s)G‘A‘"  +  W(s)G^A^)z'  + 

HI+W(s)G^)(i>-W(s)G‘AF%  O  =  x^*-A^‘x‘'  -  A  V  -  Ffx.  t). 

Proof.  A  candidate  for  the  Lyapunov’s  function 

V  =  a5(3)''3>0,  (37) 

was  built,  and  its  derivative  was  identified  on  the  states  of  the  system  (31)  as  follows: 
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(38) 


V=  3^3  =  +  W(s)G‘A“  +  W(.s)G-A^‘)z‘  +(A"  +  W(s)G‘A‘^  +  W(s)G^A-^  )z-  + 

+  (I+W(s)G-)0  -  W(s)G‘AF‘ -U+W(s)GnBhi]. 

Substituting  formulas  (35)  and  (36)  into  the  expression  (38)  the  following  expression  is  obtained: 

m 

V  =  -y-SIGN(^)-r^-^p.\3,\<0  Vp,.>0.  (39) 

Hence,  that  the  manifold  3  =  0  is  "attractive"  (or  equilibrium  point  3  =  0  is  stable),  and  the  SM 
exists  in  the  dynamic  sliding  manifold  (32),  (33)  in  the  system  (31).  ■ 

4.5  Example 

Suppose  we  have  a  nonlinear  control  system 


\x,  =  X.2 ,  Xj  =  -Xy  -x,  +Xi  fi(t),  Xj  =  X,  +  sin  X;  +  f\(t)+  lOu 
[y  =  xi 


(40) 


with  unknown  smooth  matched  I // r )  |  <5  and  unmatched  l//r)|<i  disturbances. 

We  wish  the  output  y(t)  to  asymptotically  track  a  reference  profile  y*(  t ).  That  is 
[im|y  *  (r)  -  y(0|  =  lirn|e(r)|  =  0  (41) 


The  new  basis  variables  (18)  in  the  system  (40)  have  been  introduced  as  follows: 

Zi  =  x]  -Xj,  Z2=  xl  -  x, .  Zj  =  -r*  -  .  (42) 

Assuming AF^x, r)  and  f,(t)  are  known,  the  equations  of  the  system  center  (21)  are 
written  as  follows: 

~  ^2>  ^2  ~~X2~X2+Xi  +f](t),  y  (t)  =  Xj,  (43) 

The  solution  to  the  equations  of  the  system  center  (43)  was  identified.  This  is 
xU0  =  y'(t),  x]=y*(t),  x](t)  =  -f(t)-y‘(t)-¥x]  +  f2(t).  (44) 

The  following  expression  for  the  conventional  sliding  manifold  was  obtained: 


17 


C7  =  zj-  29z.j  -  210z,  =  -'f(t)  -  SOyUt)  +  x;  -  210y'(t}  +  t\(t}  -  x,  +  29.v_,  +  210x,  .  (45) 

The  control  law  (26),  (45)  was  designed  with  (B')  'p  =  0.75 .  and  the  equivalent  control 

^eq—~^-^y  ~3y  +0.2x:,a:2 -O.l/y -2/v  +2.^i-{)Xf.,-3xj  +  l8.1x2-().\si}vCi+2.9xi.  The  system  (40) 

was  simulated  with  the  SMC  n  =  ) “‘  p  si}>n  a .  where  and  {t)  is  an 

uncertain  gain.  The  results  ot  the  simulation  show  the  perfect  output  tracking,  which  is 
insensitive  to  matched  and  unmatched  disturbances,  and  to  e  [0.7, 1.4] . 

Assuming  the  unmatched  disturbance  //r )  is  unknown,  the  simulation  was  repeated 
without  /[(f)  in  the  sliding  surface  (45).  The  output  tracking  exhibits  the  stable  output  tracking 
with  the  non-zero  output  steady  state  eiTor.  In  order  to  nullify  the  output  steady  state  error, 
assuming  Im  /,  (f)  =  const,  the  dynamic  sliding  manifold  was  designed  as 

'^  =  -y(t)-y(t)+xl-x^-W(s)e,  =  (46) 

s(s  + 50) 

The  results  ot  the  simulation  of  the  system  (40)  with  the  dynamic  SMC  (35),  (46)  show  that  the 
output  tracking  error  has  a  zero  steady  state  value  and  is  insensitive  to  the  matched  disturbance. 
/2(f).  Hence,  the  dynamic  SMC  joins  features  of  a  conventional  compensator  (accommodation 
to  unmatched  disturbances)  and  a  conventional  SMC  (insensitivity  to  matched  disturbances). 

4.6  Conclusion 

Nonlinear  output  tracking  in  MIMO  systems  with  unmatched  nonlinearities  and 
disturbances  is  considered  in  SMs.  The  system  is  nonlinearly  transformed  to  a  form  convenient 
for  SM  synthesis.  Both  conventional  and  dynamic  SMCs  are  designed.  In  both  cases,  only  the 
time-dependent  part  ot  the  unmatched  disturbance  must  be  estimated.  Results  show  that  both 
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controllers  are  insensitive  to  match  disturbances;  however,  dynamic  SMC  accommodates  to 
unmatched  disturbances.  Future  research  will  involve  the  chattering  elimination  and  discrete 
realization  of  the  designed  SMCs,  as  well.  The  developed  approach  will  be  extended  to  a 
nonminimal  phase  nonlinear  output  tracking.  Preliminaiy  results  look  promising. 

5.  Finite-reaching-time  continuous  SMC  for  MIMO  nonlinear  systems  [17] 

Sliding  mode  control  is  a  nonlinear  robust  control  technique,  which  is  actively  developing 
during  last  30-35  years  [1,2].  Since  the  closed-loop  system  exhibits  he  desired  robust’ 
performance  moving  on  the  sliding  surface,  it  is  highly  desirable  to  reach  the  sliding  surface  in  a 
finite  time.  The  system's  motion  on  the  sliding  surface  is  called  the  sliding  mode.  Historically, 
sliding  modes  were  discovered  in  variable  structure  systems,  which  feedback  stmcture  changes 
depending  on  value  of  a  state  vector  to  drive  the  system's  trajectory  onto  the  "custom-made" 
sliding  surface.  While  on  the  surface  the  variable  structure  control  yields  an  infinity  frequency 
switching  control  [1,2],  which  is  practically  unrealizable  in  many  cases.  Dynamics  of  switching 
elements  and  finite  sampling  time  of  microprocessors  prevent  from  implementation  of  infinity 
frequency  switching  control  law  and  lead  to  so  called  conti*ol  chattering  (limit  cycles  with  finite 
frequency  oscillation  in  a  closed-loop  system).  Employing  a  two-loop  structure  of  a  sliding  mode 
controller  [14,15]  it  is  impossible  to  have  a  high  frequency  switching  controller  in  the  outer  loop. 
Outer  loop  controller  produces  a  profile  that  must  be  followed  in  the  inner  loop.  Of  course,  this 
profile  must  be  continuous.  In  the  work  [17]  a  continuous  sliding  mode  controller  is  proposed  as 
a  specific  type  of  a  sliding  mode  controller,  but  not  as  implementation  of  a  discontinuous  high 
frequency  switching  control  law.  Robustness  of  the  designed  controller,  which  is  an  issue  of 
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crucial  importance,  is  addressed  as  well.  Theoretical  results  are  validated  via  computer 
simulations. 

5.1  Problem  formulation 

Consider  the  nonlinear  square  MIMO  system 
\k  =  f(x)-¥G(x)u 

where  xe  R",  ye  R""  .u  e  R"' .  Assume  that  the  functions  jfx),  h(x)  are  smooth  vector  fields, 

and  columns  ^,  (.x:)  V/  =  l.m  of  the  matrix  G{x)e  R"^”'  are  smooth  vector  fields  as  well.  Assume 

also  that  the  system  (1)  is  completely  linearizable  in  a  reasonable  domain  A:e  T,  i.e.  the  system 
(47)  does  not  have  internal  (zero)  dynamics. 

Given  in  real  time  an  output  reference  profile  y\(t)  identify  a  continuous  sliding  mode 
controller  which 

1.  provides  a  given  eigenvalues  to  the  de-coupled  output  tracking  error  perfomance  in  the 
sliding  mode, 

2.  drives  the  closed-loop  system  trajectory  to  the  sliding  surface  in  a  finite  time. 

5.2  Basis  Transformation 

Differentiating  the  outputs,  the  system  (47)  can  be  easily  transfonned  to  a  normal 
(canonical)  form  [7] 
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_  LJh^Xx) 


+  E{x)u 


Lyh„,ix) 


where 


L,AL’fK)...L,AL'fhA 


Vi  =  l,m  are  corresponding  Lie  derivatives. 


5.3  Terminal  attractor 


A  useful  result  is  formulated  in  a  following  Lemma. 
Lemma.  The  origin  of  a  differential  equation 
z  =  -0}z‘‘,  z(0)>0,  a)>0 

is  a  terminal  attractor  with  a  finite  reaching  time 

vae(0,l). 

Indeed,  integrating  (49)  we  have  obtained 
{z( ?>)'■"  =  {z( 0 ))*"''  -{\-a)cot. 

Requiring  z(t^)=Q  implies  formula  (50). 


|£:(A:)|9£{)Vxe  F;  L'}/i,.  and  LJLyh^) 


5.4  A  finite-reaching  time  continuous  sliding  mode  controller  design 

Achieving  de-coupled  output  tracking  error  performance  in  the  sliding  mode,  the 


following  sliding  surfaces  are  designed 


(52) 


(J.  +C:. 


e:  . 


Vi  =  l.m. 


(l>g^  _  _ 

where  e,  =  v/0»  e,''’  =  — — ^  and  the  coefficients  c, ,  V/  =  l.mV/  =  ()./•.  -2  are  chosen 

dt‘  ’ 

to  provide  given  eigenvalues  placement  to  the  de-coupled  differential  equations  cr,  =0 
Vi  =  l,m. 


Designing  a  continuous  sliding  mode  controller  for  the  system  (48)  we  need  to  provide 
asymptotic  stability  to  its  motion  in  the  tr  -subspace,  which  is  described  as  follows: 

(t  = '?(.)- (53) 

where 


'!'(.)  =  Vi  0,^2  V'™  (.)F .  Wi  (•)  =  y):;'  (0  + 


(I) 


-  (a:)  Vi  =  1,  /n .  (54) 


Introducing  a  new  control  variable 

u=E{x)u,  (55) 

the  system  (53)  can  be  rewritten  in  a  scalar  format.  This  is 


d.=vr,  (.)-«,  Vi  =  l,m. 


(56) 


The  system  (56)  is  obviously  input-output  de-coupled.  Hence,  analyzing  stability  of  the 
origin  of  the  system  (56),  we  can  analyze  stability  of  each  i***  equation  V/  =  l,m.  Candidates  to 
Lyapunov  functions  are  introduced.  They  are 


(57) 


We  will  look  for  the  derivative  of  the  candidates  to  Lyapunov  functions  (57)  in  the  following 
format: 

7;  =  -Q)yt‘,  a,-  >  0,  tu,>0 .  (58) 
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Then,  applying  Lemma  we  can  conclude  that  the  origin  of  the  i"’  equation  of  the  system  (56)  will 
be  reached  in  a  finite  time  V/  =  1,/n 


(v^ro)r- 


Va,  e  (0,l) . 


(59) 


It  means  that  cr^  -  0  V/  =  l,m  become  terminal  attractors.  Substituting  expressions  (52)  and  (56) 


into  formulas  (57),  a  control  function  u  is  identified  as  follows: 

«=-TO  +  /?I.  (60) 


0): 


where  R  =  diag  <{  ^  j-  V/  =  I,  m ;  Z  =  {jf"' 


CT-*  ' 


f. 


Substituting  equation  (60)  into  equations  (56),  we  obtained  equations  of  the  closed-loop 
system’s  (48)  motion  in  the  cr -subspace.  Tliey  are 

=  -2'“'  £y,.£T-“'''  V/  =  .  (61) 

Making  coordinates  of  the  control  function  (14)  continuous  and  avoiding  singularity  in  the  origin 
of  the  system  (61),  we  will  look  tor  the  parameters  "a/*  in  the  following  foimat: 


a.  =-£l  V/  =  l.m, 
(li 


(62) 


where  and  q.  are  integer  numbers.  Since  cr,  can  be  either  positive  or  negative  the  values  of 

and  2p.-q.>Q  must  be  odd.  Hence,  the  parameters  p,and  q^  must  meet  the  following 
conditions 

€  Q,:{p,.,^i|2p,>^„  p.<q.,  q.  is  odd}  V/  =  l, in.  (63) 


23 


As  soon  as  the  values  ot  the  parameters  /?,  and  q.  have  been  identified  in  accordance  with  the 
conditions  (63)  the  value  of  the  parameters  o),  must  be  chosen  to  achieve  a  given  reaching  times 
(59). 

The  original  finite-reaching-time  continuous  sliding  mode  control  law  u  has  been 
identified  in  accordance  with  foimulas  (55)  and  (60)  as  follows: 

w  =  £:(x)‘‘[-'F(.)  +  /?z].  (64) 

The  obtained  results  are  formulated  into  the  following  Tlieorem. 

Theorem.  Tlie  sliding  suitaces  (56)  will  be  reached  in  finite  times  (59)  via  the  continuous  sliding 
mode  controller  (64)  if  values  of  parameters  "a,  "  are  chosen  in  accordance  with  formulas  (60) 
and  (63). 

5.5  Robustness 

Assuming 

4^0  =  %  (•)  +  A4^(.) ,  E(x)  =  E,  (X)  +  AE(x),  (65) 

where  A'F(.)  and  A£'(Ar)  are  uncertain  vector  and  matrix  correspondingly.  Designing  control 
(64)  in  the  following  format 

u  =  E,(xr'[-%Q  +  R-'E],  (66) 

and  substituting  expression  (66)  into  equations  (53),  we  obtained 

(T  =  A'P(.)-AM-4^o(.)-(/  +  AM)/?Z,  (67) 

where  AM  =  AE(x) ■  E^^ (x) . 

Equation  (67)  is  rewritten  in  a  scalar  format.  Tliis  is 
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CO  CO 

d,=AvA,0-(AM-4^,,0),-  — .  V/  =  l./«  . 


(68) 


Substituting  equations  (68)  into  equations  (57),  we  obtained 


2-., 


£U,(1  +  AM,)  +  - 


CO  , 


-  A  V/,.  (.)  +  (AM  .  (.)),  +  X  AM,  ^<r: 


j^i 


err-' 


(69) 


Assuming 


|AM,,|<L,  <1,  |aM,|<L,,  |Ai/A,.|<i?,.  |(  AM  •%(.)),!  <G,,  |cT,  <Kj  V/ =  l,m  Vy  =  l,m  (70) 


we  introduce  positive  constants 


m  Q) 

<I>,.  =  Ri+Gi+'^—Lf7cY'~^ 

;  -l  2 
J*i 


(71) 


Transforming  equations  (69)  taking  into  account  expressions  (60)  and  (71)  we  obtained 


K<- 


2“'<l>. 


O’; 


v;."-  Vi  =  lm. 


(72) 


Inequalities  (72)  are  rewritten  in  a  format 


V,  <-S).X‘  Vi=l,m, 


(73) 


where 


^(i-«,) 


(74) 


and 
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^>  =  lfn.  (75) 

K'  I 

Inequalities  (75)  actually  form  domains  where  conditions  (73)  are  met.  It  means  that  the 
following  domains 


(T..  < 


Vi  =  1, 


m 


(76) 


a),(l-L,)-d5,  J 

will  be  reached  by  the  trajectories  cr,  in  given  finite  times  (59),  and  trajectories  cr,.  will  stay  in 
these  domains  thereafter.  Hence,  solutions  of  the  closed  loop  systems  (68)  are  uniformly 
ultimately  bounded  [8].  What  is  remarkable,  that  choosing  values  of  the  parameters  to,,  to  meet 
inequalities 


2‘''0, +d5,  — 

0);  > - ! - Vi  =  l,m 

1  La 


(77) 


and  making  2a,.  —  1  arbitrary  small,  we  can  make  a  size  of  the  attractivity  domains  (76)  arbitrary 

small  as  well,  retaining  the  control  law  (66)  continuous.  In  order  to  achieve  a  similar  result,  using 
linear/saturation  approximation  of  a  discontinuous  sliding  mode  controller,  a  gain  of  such 
controller  must  be  arbitrary  large  [7,8],  which  is  impractical.  Hence,  a  designed  continuous 
sliding  mode  controller  is  extremely  robust. 


5.6  Example 

A  nonlinear  MIMO  system  is  considered 
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X,  =  x. 


X2  =  COS-iCj  +  (.V,'  +  1)/,  +  11.  +  /,  {t) 

=  sin  +  0.5i(,  +  (.Vj  + 1)/2  +  /2(0 
Vl  =  •’^1  ’  >’2  =  -^2  . 


(78) 


Functions  /,.  (0 :  |/,.  (f)|  <  77,  V/  =  1,2  are  unknown,  bounded,  matched  disturbances.  The  system 
(78)  is  completely  linearizable  with  a  vector-relative  degree  r  =  {2,  l}.  The  system  (78)  is 
transfonned  [1,2,7]  into  the  normal  fonnat  (48)  as  follows; 


— 

COS.X3 

+ 

7.(0' 

-h 

.^2. 

sin;c. 

72(7 

(vf  -m)  1 

L  0.5  fe+l)l»2j‘ 


(79) 


The  goal  is  to  design  a  robust,  continuous  sliding  mode  controller  to  provide  de-coupled  tracking 
of  smooth  output  reference  profiles  (t)  and  y.,  (r)  in  fmite-reaching-time  sliding  modes. 

The  sliding  surfaces  are  designed  in  the  format  (52)  as  follows 

(7i=ei+5e^,  (T2=e2,  (80) 

where  e,  =  -  y, ,  ej  =  y,,  -  yj . 

The  parameters  /?,. ,  ^,.and  o),.  are  identified  Vi  =  1.2 .  They  are 

p,  =^3=2,  =(^2=3,  CO,  ==16,  CO3  =10. 

The  continuous  "cubic  root"  sliding  mode  controller  is  designed  in  the  format  (64)  as  follows: 


xl  -fl 

1 

(jcf  +lX-^2  +l)-0.5 

(x,^-l-lXx3 -hl)-0.5 

i7, 

0.5 

+1 

3. 

+1X^:3+!)- 0.5  1 

;a:,‘ +1X4+1)- 0.5  _ 

(81) 


where 


M,  =  yv  +5y,,  -cos;c3  -5X3  + 10.0^/^ 
“2  =  yir  -sinjr,  +6.3\f^ 


(82) 
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The  system  (78)  was  simulated  with  the  continuous  "cubic  root"  sliding  mode  controller  (82)  and 

with  the  continuous  "linear"  sliding  mode  controller 

/“i  =  -cos  ATj  -Sjc,  +  16(Tj  , 

1“2  =  V2r  “ sin 

It  is  worth  to  note  that  gains  before  the  cr  -  terms  in  formula  (83)  are  increased  from  10  to  16  and 
from  6.3  to  20  comparing  with  formula  (82),  which  helps  the  “linear”  controller  (83)  to  be  more 
aggressive. 

The  following  initial  conditions  and  the  reference  profiles  were  used 
x/  0 )  =  0.2,  jc,!'  0 )  =  0. 1,  JC3  (0)  =  -0. 1 5 ,  yjt)  =  0.5  sinr,  yjt)  =  0.5  cosf .  (84) 

Simulations  confirm  a  finite-reaching  time  property  of  the  designed  continuous  sliding 
mode  controller  and  show  improving  robustness  comparing  with  “linear^’  sliding  mode 
controller. 


6.  Concepts  for  a  Control  System  Development 

Engineering  analysis  is  perfoimed  to  investigate  the  implementation  of  a  sliding  mode 
controller  algorithm  that  will  reduce  sensitivity  to  feedback  noise  in  Control  Actuation  Systems 
(CAS)  such  as  those  used  on  Guided  MLRS,  Low  Cost  Precision  Kill  (LCPK),  and  Compact 
Kinetic  Energy  Missile  (CKEM).  The  GFE  model  of  the  CAS  is  used  to  demonstrate 
performance  enhancement  achieved  with  the  sliding  mode  controller  for  a  variety  of  potential 
noise  sources.  Two  control  systems  are  designed. 

First  The  “nested”  structure  of  the  control  system  is  developed.  Employing  backstepping 
concepts  [16]  the  outer  loop  controller  is  designed,  in  particular,  as  a  continuous  finite-reaching 
time  SMC,  using  continuous  sliding  mode  control  technique  [7-9,17].  Armature  current  is 
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treated  as  a  virtual  control  in  the  outer  loop  sliding  mode  controller.  Armature  current  command, 
which  is  formed  in  the  outer  loop  by  a  continuous  finite-reaching-time  SMC  [17]  in  order  to 
provide  an  asymptotic  tracking  of  the  load  position  angle,  is  tracked  in  the  inner  loop.  The  inner 
loop  controller  is  designed  in  order  to  follow  armature  current  command.  The  inner  loop  high 
frequency  switching  SMC  is  designed  in  terras  of  armature  voltage.  This  design  concept  is  based 
on  assumption  that  compensated  dynamics  of  the  inner  (current)  loop  is  much  faster  then  the 
compensated  dynamics  of  the  outer  (position)  loop.  Advantage  of  the  two  loop  SMC  is  that 
sliding  surfaces  are  of  reduced  order  and  the  disturbance  F(.)  is  matched  [1,2].  It  means  that  the 
SMC  will  compensate  for  the  disturbance  F{.)  without  its  measurement  or  estimation. 
Implementation  of  the  SMC  is  considered  in  a  relay  format  and  in  a  pulse-width-modulator 
(PWM)  format. 

Second.  Single  loop  SMC  is  designed  using  the  system  center  design  technique, 
which  is  discussed  in  the  chapter  4.  Advantage  of  this  approach  is  in  simpler  structure  (single 
loop)  of  the  sliding  mode  controller.  The  single  loop  high  frequency  switching  SMC  is  designed 
in  terms  of  armature  voltage.  Implementation  of  the  SMC  is  considered  in  a  relay  format  and  in  a 
pulse-width-modulator  (PWM)  format.  Sliding  mode  observer  is  designed  to  estimate 
disturbance  F(.)  that  is  unmatched  in  this  case. 

7.  Two  loop  sliding  mode  controller  design 
Designing  a  two-loop  SMC  for  the  servomechanism  (1),  we  start  with  the  outer  loop 
continuous  finite-reaching-time  SMC  design. 

7.1  Outer  Loop  Continuous  Finite-Reaching-Time  SMC  Design 

The  mathematical  model  describing  the  "armature  current  -  load  position"  relationship  is 
derived  from  equation  (1).  This  is 
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(85) 


(It 


dco 


=7k-,-P()l 

at  J 


y  =  e, 


Assuming  the  load  position  reference  profile  is  given  in  real  time,  the  problem  is  to 

design  a  continuous  virtual  sliding  mode  controller  (t)  achieving  an  as)miptotic  tracking 
lm|6>^rr)-6>^(0|  =  0.  (86) 

At  the  tirst  step  ot  the  outer  loop  SMC  design  the  sliding  surface  is  identified  as  follows: 
ct  =e  +  200e,  e  =  6^-9.  (87) 

While  on  the  surface  (87)  the  system  (85)  will  exhibit  the  following  tracking  behavior: 

eit)  =  eit„ )  exp(-2000 .  (88) 

that  is  insensitive  to  disturbance  F(.) . 

At  the  second  step  the  continuous  finite-reaching-time  virtual  SMC  law  is  designed  to 
asymptotically  stabilize  <t  =  0 .  The  sliding  surface  dynamics  is  described  as  follows: 

&  =9 2009^  - 1  .588cd„  +  625.2F(.)  - 278.2/,, .  (89) 

The  equivalent  control  that  achieves  the  system's  (85)  motion  in  the  sliding  surface  <T  =  0  is 
identified 

iace,  =  0.00360^  +  0.7 190^  -  0.057a)„  +  2.25F(.) .  (90) 

The  disturbance  F(.)  is  omitted  while  equivalent  control  is  estimated.  This  is 

Le,  =  0.00360^  +0.7190^  -0.057g)„  .  (91) 
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The  continuous  tinite-reaching  -time  SMC  is  designed  in  the  format  of  equation  (66)  to  achieve 
the  origin  of  the  equation  (89).  In  the  other  words,  existence  of  the  sliding  mode  will  be 
provided.  This  is 

’  (92) 

2 

where  a  =  - and  £0  =  8.  As  a  result  the  control  law  (91),  (92)  is  obtained 

=0.00360^  +0.7190^  -0.057£0,„  -5^?.  (93) 

It  is  expected  that  the  outer  loop  SMC  will  have  enhanced  robustness,  as  discussed  in  the 
chapter  5. 


7.2  Inner  Loop  SMC  Design 

The  mathematical  model  describing  the  "armature  voltage  -  armature  current"  relationship 
is  derived  from  equation  (1).  This  is 


di  1  r  1 

—  IV.  +  n,  -  RJa  -kbOiJ 

.y=ia 


(94) 


Having  the  load  armature  current  command  profile  (/)  (93)  given  in  real  time,  the  problem  is 
to  design  a  SMC  (t)  achieving  an  asymptotic  tracking 


lim|ijr)-j„(0|  =  0.  (95) 

At  the  first  step  of  the  inner  loop  SMC  design  the  sliding  surface  is  identified  as  follows: 

V=iJ0-iait).  (96) 

While  on  the  surface  (96)  the  system  (94)  will  exhibit  the  perfect  tracking  behavior: 
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Lit)  =  L(0, 


(97) 


that  is  insensitive  to  the  voltage  noise  . 

At  the  second  step  the  SMC  law  is  designed  to  asymptotically  stabilize  i’  =  0 .  The  sliding 
surface  dynamics  is  described  as  follows; 

^  =  L+  +  13.0m, „  +  277.8n,  - 277.8v, .  (98) 

The  equivalent  control  is  identified 

Vo,,  =  0.00364,  + 104  +  0.0503a), „  +  n, .  (99) 

Simplifying  the  control  law,  equivalent  control  is  estimated  at  zero  level.  The  inner  loop  SMC  is 
designed  in  the  format  of  equation  (11)  to  achieve  the  origin  of  the  equation  (98).  In  the  other 


words,  existence  of  the  sliding  mode  will  be  provided.  This  is 

v^=p  signs,  (1(X)) 

where 

p  >  |0.0036/,„  + 104  +  0.0503a)„  +  n,| .  (101) 

In  accordance  to  physical  limitation  to  the  armature  voltage  the  SMC  (100)  is  chosen  as  follows: 

=31  signs.  (102) 

Inequality  (101)  becomes 

|0.0036/^  + 104  +  0.0503a)„  +  n,|  <  31 ,  (103) 

and  gives  a  domain  of  existence  of  sliding  mode  on  the  surface  (96). 


7.3  Simulation  of  the  servomechanism  with  two  loop  SMC 
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The  system  (1)  was  simulated  with  the  two-loop  SMC  (87),  (93),  (96),  (100). 
Nonlinearities  were  taken  into  account.  Backlash  characteristic  of  the  gear  train  was  used  with 
0.75  width,  and  the  deadband  block  was  used  with  the  deadband  equal  to  3.0  (fig.  1).  The  first  set 
of  simulations  was  performed  for  =0.2iign(sint).  The  results  of  simulations  are  shown  in 

figures  3-9.  Modeling  the  noise  voltage  n,  we  assume  that  armature  current  discontinuity,  which 
is  caused  by  switching  control  voltage  v^,  leads  to  the  noise  voltage  n^..  During  simulation  the 

noise  voltage  is  modeled  as  follows:  =4-10~^-^^.  Differentiation  is  performed 

dt 

numerically. 


0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1  1.1  1.2  1.3  1.4  1.5 


Time  (sec) 

Fig.3  Armature  voltage  (0^  =0.2sign{sint)) 
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Torques :  motor  and  disturbance  (in-lbf)  Currents  (A) 


Fig.  4  Current  4  and  current  command  (6 u=  0-2sign{sin  t)) 


Fig.  5  Torques:  motor  torque  and  disturbance  torque  =  0.2sign{sint)) 
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Time  (sec) 


Analyzing  simulation  results  we  can  say  that  control  voltage  r’„  is  a  high  frequency 

switching  function  (fig.  3)  even  without  any  special  efforts  such  as  PWM.  Sliding  mode  holds  on 
time  intervals  where  sliding  surfaces  are  close  to  zero  (fig.  7).  Armature  current  accurately 
follows  armature  current  command  on  time  intervals  when  sliding  mode  holds  (fig.  4).  Load 
position  command  is  tracked  well  (fig.  6).  It  is  observed  that  due  aggressive  load  position  profiles 
(step  function)  sliding  mode  is  destroyed  on  relatively  large  time  intervals  where  the  control 
system  acts  as  open  loop.  No  effect  of  voltage  noise  (fig.  9)  to  load  position  tracking  is  observed. 

The  second  set  of  simulations  was  performed  for  =0.2sinr.  The  results  of 
simulations  are  shown  in  figures  10-16. 


4 

3.5 


■-Armature  current 
^•^ommand  to  armature  current 


3 

2.5 


•3  I - • - • - 1 - 1 - 1 - 1 - 1 _ I _ ! _ I _ I _ I _ i _ ! _ I 

0  .1  .2  ,3  .4  .5  .6  .7  .8  .9  1  1.1  1.2  1.3  1.4  1.5 

Time  (sec) 


Fig.  11  Current  4  and  current  command  4^  {&u:=  0.2sin  t ) 


Fig.  12  Torques:  motor  torque  and  disturbance  torque  (0^^  =  0.2  sin  t ) 
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Sliding  surfaces 


0  1  -2  -3  .4  .3  .6  .7  .8  .9  1  1.1  1.2  1.3  1.4  1.3 

Time  (sec) 

Fig.  13  Load  position  tracking  (9^  =  0.2 sin  t ) 
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Voluge  noise  (V)  Motor  speed  (rad/s) 


Fig.  15  Motor  speed  (0^  =  0.2 sin  t ) 
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'Ml  Mil 

■25  - 

-30 1 - ■ — I - 1 - 1 - 1 - 1 - 1 - 1 - 1 _ ^ _ I _ I _ I _ I _ 

0  1  -2  .3  .4  .3  .6  .7  .8  .9  1  l.l  1.2  1.3  1.4  1.5 

Time  (sec) 

Fig.  16  Voltage  noise  (9^  -  0.2 sin  t ) 
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Analyzing  simulation  results  we  can  say  that  again  control  voltage  v„  is  a  high  frequency 

switching  function  (fig.  10)  without  any  special  efforts  such  as  PWM.  Sliding  mode  holds  on 
time  intervals  where  sliding  surfaces  are  close  to  zero  (fig.  14).  It  is  important  that  sliding  mode 
holds  almost  all  over  the  control  interval.  Amiature  current  accurately  follows  armature  current 
command  (fig.  11).  Load  position  command  is  tracked  very  accurately  (fig.  13).  No  effect  of 
voltage  noise  (fig.  16)  to  load  position  tracking  is  observed. 

Simulations  were  repeated  with  inner  loop  SMC  implemented  in  a  PWM  format.  Dither 
signal  (pit)  was  in  a  triangle  periodical  format  with  period  equal  to  0.0001  s.  Amplitude  of  a 
dither  signal  was  equal  to  0.5.  The  inner  loop  SMC  was  implemented  as  follows: 

=31  signs,  s=s  +  ^it).  (104) 

The  results  of  the  simulations  were  pretty  much  the  same.  I  have  an  impression  that  the  voltage 
noise  does  not  significantly  affect  frequency  of  control  switching,  because  it  has  the  same 
frequency  as  the  armature  voltage  .  However,  the  switching  frequency  of  . the  armature  voltage 

Vg  looks  somewhat  higher  with  PWM  implementation  comparing  with  relay  implementation  of 

the  SMC.  It  doesn't  hurt  to  use  PWM  implementation  of  the  SMC.  And,  of  course,  it  makes  sense 
to  implement  the  inner  loop  SMC  in  a  PWM  format  to  protect  the  controller  switching  frequency 
from  other  noises.  Amplitude  of  a  dither  signal  must  be  greater  then  amplitude  of  summarized 
noise  containing  in  the  sliding  surface  equation  (96). 

8.  Single  loop  SMC  design  (method  of  system  center). 

Employing  the  method  of  system  center  (chapter  4)  we  rewrite  the  mathematical  model  of 
the  servomechanism  (1)  in  errors.  This  is 
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dQr  r*  1 

=0.00794®, .  +|a^,  -(>.00794a.,J 

dQ  ~  r 

=  350394  +  kc  -  350394,  +  78740F(.)] 

^  =  277.8(-  V,  -  lO.Or,  -  0.05035,,. )+  L  +  27784,  +  l^.Ocy,,  -  277.8n„  ] 
y^^L-^hc-^Lc] 


(105) 


where 

<^mc»  ^ac  command  signals  for  , 

~^L->  =^mc  “ »  4  =  hic  ~4  ^®  State  Variable  tracking  errors. 

Assuming  the  load  position  reference  profile  (0  is  given  in  real  time,  the  problem  is 
to  design  a  SMC  (f)  achieving  an  asymptotic  tracking 

lm|e^(r)|  =  0.  (106) 


The  command  profiles  and  £0„,  must  be  identified  through  solution  of  equations  of 
the  system  center  (22)  given  load  position  reference  (command)  profile  y,  =9^.  This  is 

=  125.961^ ,  4,  =  0.003590,,  +  2.247F(.) .  (107) 

Having  state  variable  profiles  4c  and  given  in  real  time,  the  output  tracking  problem  (105), 
(106)  can  be  transformed  to  the  following  state  variable  tracking  problem: 


d9, 

— ^  =  0.007945 
dt 

da  ~ 

— 2.  =  350394 
*  dt 

^  =  277.8(-  V,  - 10.04  -  0.0503m, „ )+  <I>(.) 
y=9. 


(108) 
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8.1  Sliding  surface  design 

At  the  first  step  of  the  outer  loop  SMC  design  the  sliding  surface  is  identified  as  follows: 

5  =  +Cid5,„ +C2e^  =0.  (110) 

In  order  to  identify  parameters  c,  and  c,  of  the  sliding  surface  (110)  the  system's  (108)  motion 
on  the  sliding  surface  (1 10)  is  derived.  Tliis  is 
do 

— ^  =  0.00794S 
dt 

=  -35039(c,<S„  +  cjj) .  (Ill) 

dt 

y=e. 

Requiring  the  desired  characteristic  equation 
(A  +  500)^  =  A“  +  lOOOA  +  250000  =  0 , 

the  corresponding  values  of  parameters  c,  and  c,  were  identified 
c,  =0.029,  C2=899. 

8.2  SMC  design 

At  the  second  step  the  SMC  law  is  designed  to  asymptotically  stabilize  5=0.  The 
sliding  surface  dynamics  are  described  as  follows: 


(112) 


(113) 


5  = -277.8V,  -17624  -6.87d5„  +«>(.). 


(114) 


Equivalent  control  that  provides  the  system's  motion  on  the  sliding  surface  is  identified  nullifying 


a  derivative  of  the  sliding  surface  (114).  This  is 

=  -6.344  - 0.025d>„,  +  ().()()360(.) .  (115) 

The  SMC  is  identified  in  the  format  (100)  as  follows: 

signs,  (116) 

where  p  >  ,  or  in  the  other  words 

p  >[-6.344  -0.0255,,,  +0.0036<1)(.)|.  (117) 

In  accordance  to  physical  limitation  to  the  aimature  voltage  the  SMC  (1 16)  is  chosen  as  follows: 

signs.  (118) 

Inequality  (117)  becomes 

j-  6.344  -  0.0255,,,  +  0.0036<D(.)|  <31  (119) 

and  gives  a  domain  of  existence  of  sliding  mode  on  the  surface  (1 10),  (1 13). 


8.3  Disturbance  observation 

Implementing  the  SMC  (110),  (113),  (118)  we  have  to  calculate  state  variable 
commanded  profiles  i^  and  m,„,,(107).  These  profiles  contain  the  disturbance  F(.),  which  can 

be  calculated  (fig.  2)  or  estimated,  in  case  when  it  is  not  known  exactly.  The  sliding  mode 
observer  [10]  was  designed  to  estimate  the  disturbance  F(.).  Designing  the  sliding  mode 
observer  we  used  the  second  equation  from  servomechanism  mathematical  model 

^  =  78740(0.4454  -  F(.)).  (120) 

Assuming  4  measuring  accessible  we  estimate  cb^  as  follows  [10]: 
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(121) 


Q),„  =-Ksign{a„,  -aj 
<Tz  =  -z-Ksign{a„,  -a)J, 

and  the  disturbance  F{.)  is  estimated  as  follows: 

F(.)  =  0.445/, -0.0000127Z. 

The  sliding  mode  observer  was  implemented  in  a  continuous  format.  This  is 

|<6,=-5OOOO(ffl.-0j 

lF(.)=0.445/,+0.635(tti. -<»,.)■ 

8.4  Simulation  of  the  servomechanism  with  a  single  loop  SMC 

The  system  (1)  was  simulated  with  the  single  loop  SMC  (110),  (1 13),  (118)  implemented 
in  a  PWM  format  (104).  Nonlinearities  were  taken  into  account.  Backlash  characteristic  of  the 
gear  train  was  used  with  0.75  width,  and  the  deadband  block  was  used  with  the  deadband  equal 
to  3.0  (fig.  1).  The  first  set  of  simulations  was  performed  for  6^^  =0.2sign{sint).  The  results  of 

simulations  are  shown  in  figures  17-23.  Modeling  the  noise  voltage  n^.  we  assume  that  armature 
current  discontinuity,  which  is  caused  by  switching  control  voltage  \\ ,  leads  to  the  noise  voltage 

n,.  During  simulation  the  noise  voltage  n  is  modeled  as  follows:  /z  =4-10"^^. 

dt 

Differentiation  is  performed  numerically. 


(122) 


(123) 
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Currents  (A)  ArmaUrre  voltage  (V) 


Fig.  18  Current  4  and  current  command  (0^^  =0.2«^n(sinr)) 
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Fig.  19  Torques:  motor  torque,  disturbance  torque  and  estimate  =  0.25i^n(sint)) 


Fig.  20  Load  position  tracking  (dj^  =  0.2«^n(sin  t)) 
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Motor  speed  (rad/s) 


-400 


.2 


Time  (sec) 


Fig.  21  Sliding  surface  O.lsign (sin  ?)) 


Time  (sec) 

Fig.  22  Motor  speed  (9^  =  Q,2sign{sin  t)) 


48 


-800 


.2 


Fig.  23  Voltage  noise  (0^^  =  0.2^i^n(sint)) 

Analyzing  simulation  results  we  can  say  that  again  control  voltage  is  a  high  frequency 
switching  function  (fig.  17).  Sliding  mode  holds  on  time  intervals  where  sliding  surfaces  are 
close  to  zero  (fig.  21).  Armature  current  accurately  follows  armature  current  command  (fig.  18). 
Load  position  command  is  tracked  very  accurately  (fig.20).  No  effect  of  voltage  noise  (fig.  23)  to 
load  position  tracking  is  observed.  The  system  was  simulated  with  the  single  loop  SMC  in  the 
relay  format  (118)  as  well.  The  results  of  the  simulations  are  pretty  much  the  same.  However,  the 
switching  frequency  of  the  armature  voltage  looks  somewhat  higher  with  PWM 

implementation  comparing  with  relay  implementation  of  the  single  loop  SMC. 

The  second  set  of  simulations  was  performed  for  0^  =0.2sinr  The  single  loop  SMC 

was  implemented  in  the  relay  format  (1 18).  The  results  of  simulations  are  shown  in  figs.  24-30. 
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0  -2  .4  .6  .8  1  1.2  1.4  1.6 

Time  (sec) 

Fig.  24  Armature  voltage  (0^  =  0.2 sin  t ) 


Fig.  25  Current  and  current  command  0-2 sin 


Load  position  (rad)  Torques  (in-lbf) 


51 


Fig.  30  Voltage  noise  =  0.2 sin  t ) 


Analyzing  simulation  results  shown  in  fig.  24  -30  we  can  say  that  again  control  voltage 
Vg  is  a  high  frequency  switching  function  (fig.  24)  without  any  special  efforts  such  as  PWM. 

Sliding  mode  holds  on  time  intervals  where  sliding  surfaces  are  close  to  zero  (tig.  28).  It  is 
important  that  sliding  mode  holds  almost  all  over  the  control  interval.  Armature  current 
accurately  follows  armature  current  command  (fig.  25).  Load  position  command  is  tracked  very 
accurately  (fig.  27).  No  effect  of  voltage  noise  (fig.  30)  to  load  position  tracking  is  observed. 

Simulations  were  repeated  with  the  single  loop  SMC  implemented  in  a  PWM  format: 
=31  sign  s,  s  =s+(p(t).  Dither  signal  (pit)  was  in  a  triangle  periodical  format  with  a  period 

equal  to  0.0001  s.  Amplitude  of  a  dither  signal  was  equal  to  0.5.  The  dither  signal  is  shown  in 
fig.  31. 
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The  results  of  the  simulations  were  pretty  much  the  same.  I  have  an  impression  that  the 
voltage  noise  does  not  significantly  affect  frequency  of  control  switching,  because  it  has  the  same 
frequency  as  the  armature  voltage  .  However,  the  switching  frequency  of  the  armature  voltage 

(fig.  32)  looks  somewhat  higher  with  the  PWM  implementation  comparing  with  the  relay 

implementation  of  the  SMC.  Also,  it  makes  sense  to  implement  the  inner  loop  SMC  in  the  PWM 
format  to  protect  the  controller  switching  frequency  from  other  noises.  Amplitude  of  a  dither 
signal  must  be  greater  then  amplitude  of  summarized  noise  containing  in  the  sliding  surface 
equation.  Implementing  the  inner  loop  SMC  in  the  PWM  format  we  could  include  an  integration 
term  in  a  sliding  surface  equation  to  compensate  for  a  steady  state  tracking  error. 
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Fig.  32  Armature  voltage  with  PWM  (0^  =  0.2sin  t ) 

9.  Discrete  implementation  of  the  SMC's 

Implementing  the  designed  SMC's  in  a  discrete  format  we  take  into  consideration  the 
following: 

1.  Command  profile  6^  (t)  has  0.0005s  computational  delay. 

2.  Censored  load  position  6^,  armature  current  and  a  motor  speed  co^  have  0.001s 
computational  delay. 

3.  Armature  voltage  (control  function)  has  0.001s  computational  delay. 

4.  The  sliding  surfaces  s  (96)  and  (110)  have  0.0001s  sampling-and-holding  computational 
delay. 
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Both  discrete  SMC's  were  simulated  in  a  relay  format  and  PWM  format.  Nonlinearities 
were  taken  into  account.  Backlash  characteristic  of  the  gear  train  was  used  with  0.75  width,  and 
the  deadband  block  was  used  with  the  deadband  equal  to  3.0  (fig.  1).  It  turned  out  that  discrete 
SMC's  implemented  in  PWM  format  showed  the  performance  with  fewer  ripples  of  armature 
currents.  The  system  (1)  was  simulated  with  the  two-loop  SMC  (87),  (93),  (96),  (100) 
implemented  in  a  PWM  format  (104)  for  =  0.2sign{sint).  A  dither  signal  (fig.  31)  was 
implemented  with  amplitude  equal  to  0.25.  The  results  of  system  simulations  are  shown  in 
figures  33-37.  The  system  (1)  was  simulated  with  the  single  loop  SMC  (110),  (113),  (118) 
implemented  in  a  PWM  format  (104)  for  6^  =Q.2sign{sint).  A  dither  signal  (fig.  31)  was 
implemented  with  amplitude  equal  to  0.5.  The  results  of  simulations  are  shown  in  figures  38-42. 
Tracking  the  load  position  command  profiles  look  good  with  both  controllers. 
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Fig.  38  Load  position  tracking  (6^  =  0.2ji^n(sin  t)\ discrete  single-loop  SMC  with  PWM 


Fig.  39  Armature  current  (6^  =  0.25/5n(sin  t)),  discrete  single-loop  SMC  with  PWM 
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tor  speed  (rad/s) 


Time  (sec) 

Fig.  40  Armature  voltage  =  0.2sign{sin  t)),  discrete  single-loop  SMC  with  PWM 


Fig.  41  Motor  speed  (9^  =  0.2sign(sin  t)),  discrete  single-loop  SMC  with  PWM 
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Fig.  42  Voltage  noise  (0^^  =  0.2sign{sm  t)),  discrete  single-loop  SMC  with  PWM 

10.  Conclusions 

Engineering  analysis  was  perfonned  to  investigate  the  implementation  of  a  sliding  mode 
controller  algorithm  that  will  reduce  sensitivity  to  feedback  noise  in  Control  Actuation  Systems 
(CAS)  such  as  those  used  on  Guided  MLRS,  Low  Cost  Precision  Kill  (LCPK),  and  Compact 
Kinetic  Energy  Missile  (CKEM).  The  GFE  model  of  the  CAS  was  used  to  demonstrate 
performance  enhancement  achieved  with  the  sliding  mode  controller  for  a  variety  of  potential 
noise  sources.  Two  types  of  sliding  mode  controllers:  two-loop  SMC  and  single  loop  SMC  were 
designed  addressing  the  servomechanism  output  tracking  problem.  Both  SMC  demonstrated  high 
accuracy  following  of  load  position  command  profiles.  Robustness  to  disturbances  and  noise 
voltage  was  demonstrated  as  well.  Control  voltage  operates  in  a  high  frequency  switching  mode 
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without  any  special  efforts  such  as  PWM.  However,  simulations  of  discrete  SMC's,  that  are 
implemented  in  a  PWM  format,  showed  fewer  ripples  of  armature  currents.  Also,  it  is  preferable 
implementing  the  designed  SMC's  in  the  PWM  format  to  protect  the  controller  switching 
frequency  from  various  noises.  Amplitude  of  a  dither  signal  must  be  greater  then  amplitude  of 
summarized  noise  containing  in  the  sliding  surface  equation.  Implementing  the  SMC's  in  the 
PWM  format  we  can  include  an  integration  term  in  a  sliding  surface  equation  to  compensate  for  a 
possible  steady  state  tracking  error. 
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